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C^ ' In the kinetic theory of dense fluids the many-particle collision bracket inte- 

gral is given in terms of a classical collision operator defined in the phase space. 
To find an algorithm to compute the collision bracket integrals, we revisit the 
eigenvalue problem of the Liouville operator and re-examine the method pre- 
viously reported [Chem. Phys. 20, 93(1977)]. Then we apply the notion and 
concept of the eigenfunctions of the Liouville operator and knowledge acquired 
^L ■ in the study of the eigenfunctions to obtain alternative forms for collision inte- 

f~^ I grals. One of the alternative forms is given in the form of time correlation func- 

CO ' tion. This form, on an additional approximation, assumes a form reminiscent of 

*^ I the Chapman-Enskog collision bracket integral for dilute gases. It indeed gives 

ly^ ■ rise to the latter in the case of two particles. The alternative forms obtained 

^P I are more readily amenable to numerical simulation methods than the collision 

bracket integras expressed in terms of a classical collision operator, which re- 
quires solution of classical Lippmann-Schwinger integral equations. This way, 
the aforementioned kinetic theory of dense fluids is made more accessible by 
numerical computation/simulation methods than before. 



I. Introduction 

Classical collision operators^^ defined in the phase space are classical me- 
chanical analogs of quantum mechanical collision operators defined in terms of 
Liouville-von Neumann operators and the classical limits of the latter. They 
appear in the statistical mechanical formulas of transport coefficients in ki- 
netic theory of matter and present an important problem to resolve in the final 
stage of implementing the dense fiuid kinetic theory to make it connect with 
experiments. Although practical methods of computing the classical collision 
operators are essential to making kinetic theory a useful molecular theory of 
matter, they have not been given much attention beyond the formal theory 
level in contrast to their quantum mechanical analogs in the literaturepEEl and 
computation of them and related expressions appearing in kinetic theory still 
poses a theoretical challenge in nonequilibrium statistical mechanics. In this 
paper, we would like to take up the subject for study and make an attempt to 
mitigate the situation. 

Classical collision operators are generally defined in formal analogy^ to 
quantum mechanical collision operators in Hilbert space, but their definitions 
have been made without much attention paid to the space of functions in which 
the operators live. Clearly, study of classical collision operators would require a 
space of eigenfunctions for the underlying Liouville operator in the phase space. 
However, the eigenfunctions of the Liouville operator are little known except 
for free particles. Here we first consider the classical eigenvalue problems for 
Liouville operators, especially when the spectrum is continuous, and then ap- 
ply the eigenfunctions to recast many-particle collision bracket integrals, which 
are given in terms of classical collision operators and appear in the kinetic the- 
ory of dense fluids. In Ref. 9, the present author made a study of eigenvalue 
problems for Liouville operators from the prespective of scattering theory of 



few particle systems. In the present paper, we revisit the problem and improve 
upon the method of calculating eigenfunctions before applying the notion and 
existence of eigenfunctions to examine a way to calculate the collision bracket 
integrals in the context of the kinetic theory of transport processes, particu- 
larly, in liquids. We will examine the eigenfunctions for classical scattering of 
particles in phase space and make use of the acquired results to cast the colli- 
sion bracket integral given in terms of a classical collision operator into a form 
more readily amenable to numerical computation/simulation methods. Since 
the many-body problems involved in such a theory do not allow simple analytic 
forms for the collision cross sections and transition probabilities for physically 
realistic dynamical events and hence pose a barrier to overcome and thereby 
make it practicable in kinetic theory in general, it is imperative to devise al- 
ternative forms of transport coefhcients expressed in terms of classical collision 
operators, so that efficient numerical algorithms can be found for them. There- 
fore, this line of study is quite relevant to implementing the classical kinetic 
theory of matter to understand transport properties of dense gases and liquids 
studied in the laboratory. 

In Section |IT1 we re-examine the eigenvalue problem for classical Liouville 
operator and obtain the eigenfunctions in terms of a generating function for 
canonical transformation in a more concise manner than in the previous work.lSl 
The associated scattering theory in phase space is considered in Section IIIII 
The results, at least the basic concept, obtained in Section [II] will be used to 
cast the classical collision operator in a numerically convenient form in Section 
IIVI The new form of the collision bracket integral will be shown to yield its 
low density limiting form in the same form as the Chapman-Enskog results for 
the transport coefhcients for dilute gases. A form of collision bracket integrals 
involving a three-body collision operator is also presented to indicate the general 



idea of the method in the case of many-particle situations. Concluding remarks 
are given in Section [V] 

II. The Eigenvalue Problem for the Classical Liouville Op- 
erator 

In the kinetic theory of matter pEEnHHill^j^Q formalism is usually formulated 
in terms of the Liouville operator, and associated classical collision operators ap- 
pear in close analogy to quantum scattering theory. Formal theory^! of quantum 
scattering has been studied in depth from the mathematical standpoint and 
the mathematical propertieaSI^ of quantum scattering operators, such as the 
Hilbert space for the operators, are well understood at present. On the other 
hand, it cannot be said the same of the classical scattering theory based on the 
classical Liouville operator, and even the mathematical nature of the function 
space for the classical Liouville operator is not well clarified at present, almost 
all of mathematical operations involving the classical Liouville operators being 
performed by simple analogy to the quantum counterparts; see, for example, 
Refs. 2, 3, 7-11. Here in this section we will consider an aspect of the prob- 
lem, limiting our study to construction of eigenfunctions, which may be given 
in terms of the Hamilton-Jacobi characteristic function.liS The eigenfunctions 
are L^ functions normalizable to unity and have a closure in phase space. If the 
Hamilton-Jacobi equation is separable, then the eigenfunctions can be given 
in terms of quadratures. Otherwise, a numerical solution method is the only 
way left to treat the problem adequately. Nevertheless, the eigenfunctions can 
be used to formulate formal classical scattering theory in analogy to quantum 
scattering theory, and the theory thus formulated can provide a mathematically 
transparent and practicable computational algorithm for the classical collision 
problems associated with transport coefficients of gases and liquids. 

The particles are assumed to interact by the pair potentials V}fe(rjfc), where 



Tjk is the relative coordinate vector between the two particles j and k. Let 
p = (pi, • • • ^pf) and r = (ri, • • • ,rf) denote momenta and positions of particles 
with / standing for the number of degrees of freedom. The Hamiltonian in the 
relative coordinate system will be denoted by H {p,r): 



j ■> j<k 

under the assumption that the interaction potential energies are pairwise addi- 
tive. The Liouville operator is then defined by the Poisson brackets times —i 
where i = \/— T: 

With —i multiplied, the Liouville operator becomes self-adjoint. The subscript 
k denotes Cartesian components of vectors p and r. With this definition of 
Liouville operator the Liouville equation can be written as 

i—^Lp{x,t), (3) 

where x is the abbreviation for the phase x = (p, r) and p is the probability 
distribution function. Since it is possible to expand p into Fourier components 
as in the integral 

/oo 
dXexp{—iXt)'il!\{x), (4) 

-oo 

we arrive at an eigenvalue probleirP of L : 

L{x)i;x{x)^XtPx{x), (5) 

where A is an eigenvalue and ipx the eigenfunction belonging to A. This eigen- 
value problem is subject to appropriate boundary conditions. It is intimately 
related to the Hamilton- Jacobi theory in classical mechanics as will be shown 
below. 



Let us consider the eigenvalue problem, Eq. ([S]), in another form, which 
appears to be more insightful. The Liouville operator L is Hermitian in the 
space of square- integrable functions ipx, which are generally complex. Therefore, 
if ipx (x) is written in the form 

V'A(a;)=A(x)exp[ir(x)], (6) 

then the eigenvalue problem can be cast in the pair of equations 

[^,r]pr-A, (7) 

[^,^]pr=0. (8) 

Eq. ([5]) indicates that the amplitude A{x) is a function of canonical invariants. 
Therefore the amplitude A{x) is a function of the basic invariants of H; for 
exapmple, a function of the Hamiltonian and the total momentum. However, 
the precise form for the function is not obvious at this point of development. 

Since the Liouville operator is invariant to canonical transformation (p, r) — ^ 
(P, Q) where P and Q are the new canonical momentum and coordinate preserv- 
ing the forms of Hamilton's canonical equations of motion, in the new canonical 
variable system the eigenvalue problem takes the form 

L{X)MX) = \MX), (9) 

where X — (P, Q) and 

LiX) = -^[H,]p^=-iJ2( 

I. V 



dH d dH d 
^ .dPkdQk ~dQ~kdPk 



(10) 



Here H is the new Hamiltonian H — H(P,Q). If the new Hamiltonian is 
independent of Q, then we have 

fifj r) 
UX)--r[H.\,,^-,j:-—^ (II) 



We will return to this form for L presently. Other mathematical properties of 
L and aspects of the eigenvalue problem in hand are discussed in Appendix A. 
They are collected in Appendix A for completeness and also to make this article 
self-contained. 

Let us assume 5(P,r,i) to be the generating functional of canonical trans- 
formation such that 

dS ^ dS 
Pk = ^— , Qk = T7^- (12) 

OTk dPk 

The generating function S obeys the Hamilton- Jacobi equation 

Under this canonical transformation the amplitude function is given by 

A{x) =6(^(x-^y (14) 

where ai,- ■■ , a/ are the values of Pi, • • ■ ,P/, respectively, with / denoting 
the degrees of freedom. It is sufficient to take this form because the Hamilton- 
Jacobi equation is subject to a constant energy for a conservative system. It is 
easily shown that the expression for A given above satisfies Eq. ([5]) . 

To find the phase function F, we observe that Eq. ([7]), on canonical trans- 
formation, can be written as 

[iJ,r]pQ = A. (15) 

Explicitly written out, Eq. ((T5|) has the form 

^ V- 9H dT 



It is useful to define 



--!!■ <i^' 



which is independent of Q since H — H (P) is independent of Q. Then, Eq. 
P6|) is readily solved, and we find 

r^V— Qfc + ro, (18) 

where Fq is a constant, which can be absorbed into the normalization factor, 
and constants (eigenvalues) A^ are such that 



A = E^ 



fe- 

k 

They will be more precisely determined on imposing suitable boundary condi- 
tions on the eigenfunctions -ijjx. It should be noted that Qk = Qk (r, a). 

In summary of the results up to this point, under the canonical transforma- 
tion (J12p the eigenfunctions are in the form 

Afe ds 



V-A^AojJ^a-— jexpfz^- 



Ao5(a~p)expUV— Qfc), (19) 



(*?H 



where Aq is the normalization factor, and pk and Qk are used for the second 
equality of Eq. (|19p . Since the generating function S can be written for a 
conservative system in terms of the characteristic function W defined by 

S = W -Et, (20) 

where E is a, constant (energy), the eigenfunction can be also written as 

i,x = A,5{(x-^)e^v\^Y.—^\ (21) 

y ^ uik oak J 

This result shows that the eigenfunction will be explicitly found if the Hamilton- 

Jacobi equation is solved, given the initial conditions for (p, r). 

This function ipx — i/) (r, p;Q;, A) may be regarded as an eigenfunction with 

two sets of eigenvalues (a, A). This set of eigenfunctions can be easily shown 



orthogonal and normalizable. For this proof we simply note that under the 
canonical trasformation and subject to J (a — p) 

djQi,--- ,Qf) ^ d{pi,---Pf) ^ ^ 
d{ri,--- ,rf) a(ai,--- ,a/) 

Before proceeding further, it is interesting to note the following. Recalling 
that the generating function S can be expressed in terms of the action integral 



S{t)^~ f dsCip,r,s), 



(22) 



where >C(p,r, s) is the Lagrangian, the eigenfunctions, e.g., Eq. (??), may be 
expressed in the form 

(27 



tp (r, p;a. A) = - — ^6 { ct - — dsC (p, r, s) ) x 
iV^-— / ds/:(p,r,s) 



(23) 



which reminds us of the path integral"'^'* in quantum mechanics in that it is given 
in terms of the action integral. It is useful to note 

/ dsC (p, r, s) = ^ / dsp,^ ^E{t-to)^Y. ^'^^^^ - E {t - h) . 

Jto J Jto "■S j Jr,o 

(24) 
We now consider periodic boundary conditions on the eigenfunctions. Let 
n be the dimension of the cubic box. Then the eigenfunctions must obey the 
condition 

V'A(r + 0)-VA(r). (25) 

It follows 

(26) 



'^[Qj{r + ^)-Qj{r)]^2^lj. 



Here Ij is an integer. Since Qj (r + fi) — Qj (r) = 57 it follows that 



A,- 



2TrljU 

n 



(27) 



Define the wave number kj 

2ttI 



% = -^ a, -0,±1,±2,...). (28) 



Hence the eigenfunctions are: 



Ao 



tp (r, p;a, A) = - — 2_J (q, - p) exp i ^ kjQj . (29) 



i2n) 



J 



Since the characteristic function W is a, surface in the phase space and the 
trajectories of Ij ^ are on the family of surfaces, the eigenfunctions represent 
waves propagating perpendicularly to the surface of the characteristic function 
W with the phase 

3 j j 

Eqs. (|2T|) and (pO| represent the formal solution of the eigenvalue problem 
in terms of the Hamilton-Jacobi characteristic function W. This method of 
constructing the eigenfunctions is not only slightly different from, but also more 
insightful than that obtained in the previous work^ by the present author. 

For noninteracting particles we easily obtain 

V'A = AqS (a - p) exp i ^ k^Vj , (31) 

in agreement with the known result in the literature P 

It is an interesting exercise to construct the eigenfunctions for separable 
bound-state problems in the case of two interacting particles. Such examples are 
given for a few cases in Ref. 9. Such eigenfunctions can be used for calculations 
of few-body dynamics problems related to bound states or scattering of isolated 
two-particle systems in a classical approximation to quantum dynamics. 
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III. Classical Scattering Theory in Phase Space 

Collision bracket integrals appear in the formulaa^E^E^ of transport processes 
when dynamic processes are treated for gases and liquids in kinetic theory. 
If we wish to achieve a molecular theory of such processes it is essential to 
calculate them in terms of molecular information on the basis of mechanical 
laws. Since it wouldn't be possible to expect to evaluate them in analytic 
form for a realistic potential model, it would be essential to develop numerical 
computational algorithms for them. 

As a preparation for expressing collision bracket integrals involving classical 
collision operators in an equivalent but computationally more practical form 
and also to make this work self-contained, we briefly review classical scattering 
theory in the phase space. Formal classical theory of scatteringiEIMl can be 
formulated in the phase space in a manner parallel to the quatum mechanical 
scattering theory, and it holds some advantages for statistical mechanics and, 
especially, for kinetic theory as has been frequently demonstrated in kinetic 
theory investigation^^isiiilil) jji which classical collision operators are used in a 
formalism analogous to quantum scattering theory. Nevertheless, the meanings 
of such collision operators have not been studied beyond the formal definition 
level. Since they are often used in such investigations but wihtout their com- 
putational methods sufficiently well discussed, it is all the more important to 
try to comprehend their mathematical basis, so that one can perform computa- 
tions, for example, of the collision integrals appearing in the theory of transport 
processes. Our aim here, however, is not in developing classical scattering the- 
ory in depth, but rather in exposing the essential features that may be relevant 
to mathematical treatments of and developing computational algorithms and 
suitable approximations for them which one might use in future calculations of 
transport properties. For this purpose the eigenfunctions presented in the pre- 
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vious section provide valuable mathematical tools and a conceptual framework 
as well as a methodology for such efforts. 

Consider a scattering situation where particles (beams) at infinite separation 
converge toward each other and interact and then separate to infinite relative 
distance from each other. In classical mechanics, there is no concept of waves 
for particle motion. However, when the collision problem is formulated in the 
phase space, there appears a notion of waves as we have seen in the previous 
section where the eigenfunctions are calculated for continuous spectra of the 
Liouville operator. The Liouville equation 1^ governs this scattering process of 
waves in the phase space. We describe the basic aspects of the theory here. 

The collision process of the particles is assumed governed by the Liouville 
equation (jH]) over the course of collision. Since at remote past and distant future 
where particles do not interact with each other, the probability distribution 
function p also obeys the free particle Liouville equation 

*§^=LoP, (32) 

where Lq is the free Liouville operator defined by 

Pk d 
m drk 



^"— E^T^- P3) 



k 

If the system is prepared at an eigenstate A at remote past and thus we set 

p(x,t)=exp(-^At)$A(x), (34) 

the corresponding eigenvalue problem is 

Lo$A (x) = A$A (x) . (35) 

This eigenfunction is of the form as in Eq. (I5T]) . To describe the scattering 
process it is useful to introduce the interaction representation. Therefore, we 
define 

p{t) = exp {{Lot) p (t) . (36) 

12 



This puts the Liouville equation ([3]) in the form 

i^=U{t)pit), (37) 

where 

Li(f) = cxp(zLo<)Li cxp(— iLot), (38) 

with Li denoting the interaction Liouville operator 

^ drk dpk 
The formal solution for Eq. (1571) can be written in the form 

p{t) — exp(zLot) exp [— iL (i — to)] exp(— iLo<)/5(io)- (40) 

As to — ^ — oo, where the interaction between the particles vanishes, there holds 
the limit L — > Lq and therefore 

p{t) -> exp(iLot) exp(-iLoi)$A = *a, (41) 

and similarly as i — > oo. Consequently, in the limit of fo ~^ —oo the distribution 
function is given by 

p{t) = lim exp [~iL (t - to)] exp(-iLoto)/5(to) 

= lim exp[-iL(t- to)]exp(-iLoto)*J>A- (42) 

in— > — oo 

This is a strong limit which may be written in an equivalent form by using the 
Abel-Tauber theorenPEsl 

f° 

pit) = lim el ds e*^* exp [— zL (t — s)] exp(— i_Los)$A 

/■" 
= exp (— iLt) lim e / dse*^* exp (iLs) exp(— iAs)$A- (43) 

Physically, this form implies a time average of a wave train released from remote 
past until t = 0. We define a time-independent function -0^ by the formula 

W\ — lim e / dse^exp (iLs)exp(-iAs)$A (44) 
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to formulate time-independent scattering theory, since it contains all the neces- 
sary information on the scattering system. The nature of the scattered function 
tp\ can be better understood by recasting Eq. (|44l) on performing integration: 
with the definition of a complex number 

z ~ \ -\- it 

we find 

^^+^ = -ie (L - zY^ $A = -i^n (z) $a, (45) 

where TZ (z) is the resolvent operator and the limit sign e — > 0+ is omitted for 
notational brevity. The limit must be taken when the calculation of the right- 
hand side is completed. This limit will be understood henceforth. Multiplying 
(L — z) to Eq. (PSJ) from left and taking the limit e — >■ 0+, we obtain 

LV'i+^=A^i+\ (46) 

This equation suggests that ijj)^ is an eigenfunction of L belonging to the same 
eigenvalue A as for Lo . If it is recalled that the free state of the scattering system 
prepared in the remote past also has the eigenvalue A with eigenfunction $a, 
we see that the full scattering eigenfunction 'i}j\^ and the free eigenfunction $a 
belong to the same eigenvalue spectrum, that is, the scattering process occurs 
on the shell of a given eigenvalue A of the Liouville operators and the total 
eigenvalue of the system does not change even if there is scattering of particles 
by each other. This is an important point worth keeping in mind when we apply 
the classical scattering theory formalism to kinetic theory and, in particular, 
to calculation of collision bracket integrals given in terms of classical collision 
operators. 

The classical resolvent operator in Eq. (P5|) can be recast into another al- 
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ternative form 

n{z) = (L-zy^ 

= 7^o (z) - 7^o (z) Ll7^ (z) 

= 7^o (z) - 7^ (z) Ll7^o (z) , (47) 

where the free resolvent operator TZq (z) is defined by 

7^o(^)-(Lo-z)-^ (48) 

Upon use of Eq. (|T7|) . Eq. (P5|) can be written in a more familiar form 

^i+)-a>A-7eo(z)LiV'W, (49) 

which is called the classical Lippmann-Schwinger equation for the scattered 
eigenfunction ip^ in analogy to the quantum mechanical Lippmann-Schwinger 
equation for scattering.^ ^ It is an integral equation for -0^ . By defining the 
classical collision operator T (z) 

T(z)=Ll-Ll7^o(z)T(z), (50) 

we may put the classical Lippmann-Schwinger equation in an alternative form 

1^i+^ = $A-7^o(z)T(z)$A. (51) 

The operator relations (|47|) and ([50]) can be rearranged to obtain the relations 

7^ (z) Li = 7^o (z) T (z) , (52) 

Ll7^ (z) = T (z) 7^o (z) , (53) 

n (z) = 7^o (z) - 7^o (z) T (z) 7^o (z) , (54) 

T(z) = Li-Li7e(z)Li. (55) 

These relations will be found useful in various calculations involving collision 
operators since the classical collision operators appear in the theory of transport 

15 



processes. They have been in fact used in some of formal kinetic theories^EiEll 
in the hterature in 1960s and 1970s. 

The equations presented remain structuraUy the same as for both two- 
particle and many-particle collisions, and since they can be easily transcribed 
into many-particle versions, we will simply consider the equations presented here 
applicable to a many-particle situation. Since the many-particle collision oper- 
ators are often computed in terms of collision operators of a smaller number of 
particles — usually two-particle operators — owing to the inherent difficulty with 
many-particle dynamical problems by using some sort of expansion, there are 
some precautions that must be taken when many-particle collision and resol- 
vent operators are calculated in terms of operators involving smaller numbers 
of particles, for example, two-particle, three-particle operators, and so on. For 
example, the often used binary collision expansiorPEiEIl fQj- many-particle col- 
lision operators is accompanied by a divergence difiiculty, which can be avoided 
by using a suitable cluster expansionj ^ l ^^ l Since many-particle operators are of- 
ten decomposed into such lower-order operators of smaller numbers of particles, 
the precautions mentioned become important when kinetic theory results are 
computed for experimental comparison. Detailed discussions on these points 
are given in Chapter 9 of Ref. 3 to which the interested reader is referred. 

IV. Alternative Forms for the Collision Bracket Integral 

In kinetic theory, either the kinetic equation is derived approximately from 
the Liouville equation in a form breaking the time-reversal invariance or an ir- 
reversible equation is postulated on the basis of the viewpoint that the kinetic 
equation is a fundamental and irreversible equation for mesoscopic description 
of a many-particle system, built on (reversible) classical or quantum mechan- 
ics, not something that canbe derived from time-reversal-invariant Hamiltonian 
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equations of motion. The difference between these approaches is just philo- 
sophicaL Either way, we are treating macroscopic systems from the molecular 
viewpoint by using an irreversible equation, which is qualitatively different from 
the Newtonian equation of motion or the Schrodinger equation in the sense that 
the former is time-reversal symmetry breaking whereas the latter is time-reversal 
invariant. The time reversal symmetry is broken by the kinetic equation (e.g., 
Boltzmann equation) because of the time-reversal symmetry-breaking collision 
term (e.g., the Boltzmann collision term) in which energy dissipation accompa- 
nying irreversible processes is vested. The collision term is given in terms of the 
collision cross section or collision operator which is relatively simple to calcu- 
late in the case of two-particle elastic scattering, but becomes hard to calculate 
as the complexity of the system increases if the number of particles involved 
increases beyond two. All material functions such as transport coefficients are 
given as collision- weighted averages of dynamical quantities, and the theoretical 
formulas for material functions derived by any kinetic theory cannot transcend 
the formalistic level to make connection with experimental data unless a practi- 
cable method of computing collision- weighted averages is developed. By method 
we do not mean some sort of analytical theoretic method, but an algorithm by 
which the aforementioned averages can be efficiently computed on a computer, 
for example, in a manner similar to the molecular dynamic or Monte Carlo sim- 
ulations for equilibrium fluids. Therefore, the aim is to transform the formal 
expressions for material functions in kinetic theory to relatively simple quadra- 
tures or forms readily computable electronically. It is the main objective in this 
section and, in fact, of this paper. 

To the end stated, it is sufficient to consider the following collision bracket 
integral appearing in the dense fluid kinetic theory formulated elsewhere (Ref. 
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3): 

Ic = i {h{x)T (z) h{x)Feq (x)) = iJc, (56) 

where h{x) is a function of phase x — (p, r) = (p'^^jr*^^^) for an N particle 
system; for example, in the case of dilute gases 

h{x) = Y. "^^- [C.C,]^') , \ Y. "^^-^'C, (C, = V, - u) (57) 

3 3 

with Vj and u denoting respectively the particle velocity and the mean fluid 
velocity; Feq{x) is the equilibrium distribution function; T (z) is the iV-particle 
collision operator obeying the iV-particle version of the classical Lippmann- 
Schwinger equation (1501) : and i = a/— T. The first member on the right hand 
side of Eq. ((57)) is for the shear stress and the second member is for the heat 
flux in the case of a dilute monatomic gas. In the case of dense fluids, h{x) 
consists of virial tcntors, energy flux tvectors, etc.; see Ref. 3 for explicit forms 
for them and also Eq. (|82|) for the virial tensor below. The symbol [B] stands 
for the traceless symmetric part of second rank tensor B. The angular brackets 
stand for integration over the phase space. The limit e — > 0+ must be taken on 
completion of calculation for the average in Eq. ([5S)) . We aim to cast the phase 
space average in Eq. (j56p in an alternative, and more readily computable, form 
since it is not clear at present how the classical collision operator T(z) may be 
numerically evaluated on a computer. 

It is convenient to use the bra and ket vector notation so that we can write 

Je = {9\T (z) I g) , (58) 

where 

\g) = \h{x)F^/^). (59) 

Let <i>A denote the complete set of eigenfunctions for Lq, the free Liouville op- 
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erator for the N particles. The complete set has the closure 

J2\^,ix)){^,{x)\=l. (60) 

A 

The collision integral J^ then can be written as 



Jc^J2J2(9\^x(.^)) {^x(.^)\T{z)\^X' (x)) {<^y ix)\g) . (61) 

A A' 

Upon rearranging Eq. (|54p into the following form 

T(z) = 7^o-l (z) [7^o (z) - n (z)] n^^ (z) (62) 

and using it in Eq. (|6ip. the collision integral can be written as 

•^^ - E E (^ - 'T (^' - ^) (ff l*A (x)) ($, (x)! n (z) |$v (2;)) ($v ix)\g) 

A A' 

- E E (^' - ^) (-9 I'^A (^)) (*> (^) I '^A' (x)) ($A' (x) I 5) . (63) 

A A' 

Since the eigenfunctions are orthogonal, the collision integral is further simplified 
to the form 



^^ = E E (^ - ^)* (^' - ^) (-9 l*A (x)) (<J>A ix)\ n (z) |$v (x)) (<J>v (x)| .9) 

A A' 

^J2{X~z){g\<^,{x)){<^,{x)\g). (64) 

A 

The scattering theory consideration made earlier in the previous section indi- 
cates that the collision occurs on the shell of a given A. This means that there 
will be negligible contribution from the off-shell A and thus the important con- 
tribution is made from the term where 

{X - z)* {X' - z) = - {tef = e\ 

Therefore, we obtain 

^= - ^' E E (-9 1'5'A (^)> i'^x (^)l 7^ (z) l*A' (:.)> (^A' ix)\ g) 

A A' 

-Z6^(g|$;,(x))($A(x)|ff) 
A 

= €^g\'R{z)\g)~i€{g\g). (65) 
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By using the definition of the resolvent operator the integral Jc can be written 
in the form of integral 

/"OO 

.h = ie dte-''{e'^'g\g)-ie{g\g). (66) 

Therefore, the collision integral Ic may be written as 



Since 



{g\g}-e dte'^' (e^'^'g \ g) 
Jo 



git) = e^^'g = F,'/'e^% 



(67) 



owing to the fact that the equilibrium distribution function consists of invariants 
of L, we finally obtain 

Ir = e {g\g)-e I dt e-^* {g{t) \ g) 
Jo 



dte-"' 



{git) - 5 I ff) 



(68) 



This form and, particularly, the second term on the right in the first line, is 
rather reminiscent of the time autocorrelation functions appearing in the linear 
response theory PSIIlEIl fjere e may be taken with the inverse of a sufficiently 
large value of time, Tc, so that 



where 



dte-*/^^{git)-g\g) 
dse-^Ag{sT,)\g) {s = t/T,). 



^g{sTc) ^ g- g{sTc). 



(69) 



(70) 



This is an alternative form for the collision bracket integral in Eq. (I55|) . It is 
given in terms of a time correlation function which is certainly more readily 
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amenable to molecular dynamic simulation methods than the one involving the 
classical A^-particle collision operator T(z). 

If there exists a plateau value region in the integrand, lS.g (sTc) — g — g{sTc) 
may be approximated by /S.g (tc), the integral may be approximated by 

/, = l(Ag(r,)|g), (71) 

Tc 

By reversing the collision process in the phase space and using the time reversal 
invariance of the phase volume, it is possible to put this integral in a symmetric 
form 

/, = -L(A.g(r,)|A3(r,)) 

ZTc 

= ^J d.(^) A.g (r,) A.9 (r,) F^^ (.(^)) , (72) 

which is reminiscent of the collision bracket integrals in the Chapman-Enskog 
theorjPSED gf dilute gases, but this form is not limited to dilute gases. This is 
another alternative form for the collision bracket integral we set out to show in 
this paper. 

This integral indeed can be shown to be the Chapman-Enskog collision 
bracket integral for Ag (tc) if we consider the case of two particles for Ic- For 
this purpose, let us write I^ explicitly for two particles: 

1- J d.(2) A5 (r,) Ag (r,) F^^ (.t^)) . (73) 

If it is assumed that Ag {tc) depends only on the momenta and Feq {x^ ) = 
feq{Pi)feq{P2) as is the casc of the Chapman-Enskog collision bracket integrals, 
then it is possible to show^^'^'^^ that 

dx(2) ■■■ = TcV dvi / rfv2 / dip dbbgi2 ■ ■ ■ (74) 

in the standard notation. Here V is the volume of the container, b is the impact 
parameter, ip is the azimuthal scattering angle, gi2 is the relative speed, and 
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Ic 

2tc 



feqiPi) is a singlet equilibrium momentum distribution function. By using this 
result, we finally get 

4'^ = 2^7 ^'^i y ^^2 y d^J dbbg,2Ag (rj A5 (rj /e,(pi)/e,(P2), 

(75) 
which, apart from the constant factor, is clearly the Chapman-Enskog collision 
bracket integraPSMI fgj. /\^g (^-^^^ = g — g {tc) where g (tc) is the post-collision 
value of g if t^ is taken for the collision time for the particles. 

In the case of Ic involving three particles, we employ, for example, mass- 
normalized coordinates^ which are subsequently expressed in hyperspherical 
coordinates-l^Sl In the aforementioned mass-normalized hyperspherical coordi- 
nates the three-particle phase integral can be written apI^Zl 

/ dx^^'> = TcV I dPc I d-Ki f d-K2 I dQi f db b* {P/fi) , (76) 

where Pc is the center of mass momentum, tti and 772 are two components 
of mass-normalized relative momenta of the three particles, P = [tti +772], 
fj,^ = mira2m'i/ {mi + TO2 + TO3) = Tn^/3, 

dfl4 = sin^ 6*4 sin^ 63 sin e2deide2d93,d0i (77) 

with 6*2, 6*3, Oi denoting the hyperpolar angles (0 < 6*^ < tt : i = 2, 3, 4), Oi is the 
hyperazimuthal angle {Q < 9i < 2tt), and b is the generalized impact parameter. 
For details of the coordinate transformation and the collision dynamics leading 
to Eqs. (ffS)) and ([77]) . see Sec. 13.6 of Ref..'^ By using this result, the three- 
particle collision bracket integral can be written as 

/? ) ^l^ J d'^c J diTi J dn2 J dn^ p db b^ (P/^i) Ag (r,) Ag (r,) 

X /«g(7ri)/e,(7r2)/,,(Pc), (78) 

in the case where Ag [tc) depends on the momenta of three particles only and 

F^? = /e,(7ri)/e,(7r2)/e,(Pc). 
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(3) 

It must be noted that parameter Tc is canceled out in Ic if t^ is taken as the 
colhsion time of the three particles. Here we remark that in the case of N- 
particle collisions Ic in Eq. (j73p . that is, T (z), may be expanded into a cluster 
expansion^ where Ic and Ic appear as the leading order contributions in a 
density series. Since such an expansion is not the aim of this work, we will not 
discuss it in this paper. 

V. Concluding Remarks 

In this paper, we have re-examined the eigenvalue problem of the Liouville 
operator in the context of kinetic theory. The eigenfunctions of the Liouville 
operator are shown possible to construct by making use of canonical transfor- 
mation. The eigenfunctions obtained enable us to consider classical scattering 
theory in a mathematically precise manner. An important point realized in this 
investigation is that the classical scattering involves waves moving perpendicu- 
larly to the characteristic function W with phases given by Eq. (|30|) and occurs 
on the shell of a fixed eigenvalue which is shared by both the incident and scat- 
tered waves in the phase space. By using this latter notion in the calculation of 
the collision bracket integrals appearing in kinetic theory of dense fluids, it is 
possible to recast them in terms of a time correlation function. This latter form 
can be put in an approximate form reminiscent of the collision bracket integral 
in the Chapman-Enskog theory of dilute gases, if the there exists a plateau 
value region in the dynamic quanties involved in the time correlation functions. 
Otherwise, the collision bracket integrals are expressed in the form of a time 
integral of time correlation functions of Ag (sTc) averaged over the equilibrium 
ensemble distribution functions. 

The recast forms for the collision bracket integral are certainly more suitable 
for numerical computation on a computer than the form given in terms of the 
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classical collision operator T (z). The two-particle example is shown to give rise 
to the traditional Chapman-Enskog theory result for dilute monatomic gases. 
We have also presented a three-particle collision bracket integral in the case of 
dilute gases where there is no statistical correlation. The alternative forms of 
the collision bracket integral present a possibility of developing a new way of 
computing many-particle collision bracket integrals on a computer. With this 
investigation we now have achieved a numerical algorithm to compute the trans- 
port coefficients appearing the kinetic theory of dense fluids that is comparable 
to that of linear response theory. 

Finally, we note here the relation of the collision bracket integral to the 
transport coefficient. For example, the viscosity is related to the collision bracket 
integrals in the following manneiC^ 

5m 1 

where a is the diameter of the molecule, m is the reduced mass, n is the number 
density, and 

r ^ 



j,k 



.-1 



T^ 



ds{Ah{sTc)\Ah{sTc)), (80) 



Ah (sT,) - Y. H'^ (^^-) - ^r (0)] (81) 



with h- standing for the virial tensor 



duji{rji) 



ft^.^^ = [rnC,C,f^^-±^/[r,,r,,]^^^ (82) 



ji - — ^ rji = |rj - r;| , Uji = potential energy (83) 

arj 

1.. ... 1 



rji 



[A]^'^' = - (A + A*) - -^^TrA, traceless symmetric part of tensor A. (84) 

which is in the form, namely, the collision bracket integral considered in the 
previous section. The collision bracket integral on the right can be expressed in 
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the alternative form, Eq. (|72p . presented in the previous section. Thus the vis- 
cosity is inversely proportional to the collision bracket integral. Separating the 
center of mass part of dx^^'> so that dx^^^ = drdPdx^^/''\ {Ah (sTc) \Ah (sTc)) 
may be cast into the form 

{Ah (sT,) \Ah (sT,)) ^ Jdpf da;(^/=)Fif/^) (a;(^/^)) Ah (sr,) : Ah (st,) , 

(85) 
where Feq is the equilibrium distribution function excluding the center of 
mass motion part of the N particles. 

Other transport coefhcients may be expressed similarly in terms of the col- 
lision bracket integrals considered earlier in this work. It should be noted that 
the manner in which the time correlation function Eq. (|72p appears in the 
transport coefficient and also the meaning of the time correlation function are 
different from those in the linear response theory. This difference stems from 
the fact that the present kinetic theory and the linear response theory differ in 
their basic approach to transport processes. 
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Appendix A. Some Mathematical Properties of Eigenfunc- 
tions and Eigenvalues 

The Liouville operator L may have eigenfunctions and eigenvalues, but also 
has some important mathematical properties, which are useful and necessary 
for developing physical theory with it. We list them in the following. 

Property 1. The Liouville operator is hermitean. Consider 

^AA' = / dx'ipl,Lipx. 

Then 

lly = / dx (tPyLi^x)* 

= I dxipx'Li* 

dx-4)yi[H,-)pl]j^g 

= I dxrx{-i)[H,M^^ 

dx^/jl'Ltpx' 
~ lyx 
Therefore L is hermitean and hence 

V = L. (Al) 

Propertiy 2. Since L is hermitean, its eigenvalues are real. Let 

Lipx = Xipx- 
Then 

L*rx = >^*rx- 

Hence it follows that 

dx irx^^Px - ^xl^*rx) - (A - A*) / dx^x^x 
= A-A* 
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But by virtue of hermiticity 

/ dx (V'aLi/'a - V'aL*V'a) = 0- 

Therefore 

A = A*, (A2) 

which imphes A is real. 

Property 3. The cigcnfunctions belonging to different eigenvalues are orthog- 
onal. Proof for this follows from Property 2. 

Property 4. Any function in the phase space can be expanded into the eigen- 
functions of L. 

f{p,<i)^^aktpkip,(l), (A3) 

k 

where 

ftfc = / dTf (p, q) tpl (p, q) . 

This property follows from the fact that L is a linear operator. In the case 
of time-dependent functions 

p(p,q,t) =^afc(i)Vfe(p,q). (A4) 

k 

Inserting it into the Liouville equation yields 



.dp 



k 

^ flfe (i) AfefAfc (p, q) 

k 

dak 
~dt 



^^-7^■^k{p,^), 



k 

which gives rise to the ordinary differential equation 



i—— = Afeflfc (i) . (A5) 

dt 
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Integrating this wc obtain 

Ofe (t) = Ck cxp {-iXkt) 

Finally, 

p(p,q,t)=^Cfee-^^'=Vfe(P,q)- (A6) 

fe 

Of course, the eigenvalues and eigenfunctions appearing here are not 
known in terms of the system properties. Nevertheless, it is possible to 
deduce some general properties of the expansion, eigenvalues, and eigen- 
functions, and they may help assess eigenvalues and eigenfunctions we will 
obtain later. 

Property 5. Since the distribution function is normalized, 

J dxp{p,q,t)^l, (A7) 



it follows 



Since 



k 

j dxL,-)pk (p, q) = Afc / dxipk (p, q) 



but 



/ dxLtpk (p,q) = -i'Yl dx 



fdHd^k dH dij) 



J „ \dpj dqj dqj dp J 



i"^ dx [ TT— HI c'^/a.. ] ^k 



J 



OH OH 

•iidpj dpjdqj 



' / ^ / ^^ I Q y^k\qj— boundary ^^ Yklpj— boundary 



if 







o Wk\qj—bou7idary — ^ Wk\pj— boundary U, 
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it follows 

\k [ dxi;k{p,q) = 0- (A8) 

This means 

/ dxiJk (p, q) = (A9) 

if \k 7^ 0. However, A^ 7^ is not equal to zero for /c = 0. Therefore the 
integral of eigenfunctions for fc 7^ vanish. It thus follows 



Co / dxxpo{p,q) == 1. 



That is, 

co = l/ Jdxi,oiP,ci) (AlO) 

and the normalization of p (p, q,t) is fulfilled. 

In this manner, the distribution function in the nuclear space — the phase 
space — is expandable in the eigenfunctions and its expansion is reconciled 
with the notion of eigenfunctions being C^ functions in the Hilbert space. 

Property 6 It is well known that the Liouville operator is invariant to time 
reversal times complex complex conjugation. That is, if we denote by 9 
the time reversal operator then 

6IL = - L. 

Therefore if we denote the compound operator of time reversal and com- 
plex conjugation by 'd 

'& = 9{*) = {*)e 

then 

Therefore the time reversed distribution function p (— p, q, — t) obeys the 
same Liouville equation as p (p, q,i). We note that this is in the root cause 
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of the difficulty in kinetic theory of matter when we develop it by means 
of the Liouville equation, which is time reversal invariant. In connection 
with this operator •&, we note the following: Since 



and the eigenvalues are real we have 

L>I. = AfeV^I^ (All) 

and hence 

This means that f/;^ is the eigenfunction belonging to —Xk- Furthermore, 
since 

6*2 = 1, 

it also follows 

OL* = -61, = i?L = L. 

Now operating 9 on Eq. (JA11[) we obtain 

but 

eL*rx, = L0^*^ 

and hence 

Lorx, = Xkorx,, 

which implies that 

Oipl^ = an eigenfunction of L belonging to A^ = ipx^, . 

That is, 

V4-=^V'A.. (A12) 

This property is important for proving the orthogonality of eigenf unctions. 
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Property 7. The eigenfunctions '0a are generally complex and hence they can 
be written in terms of amplitude and phase: 

t/^A (p, q) = Ax (p, q) exp [iVx (p, q)] (A13) 

where 

LAA(p,q)-0, (AM) 

LrA(p,q) = -a. (A15) 

This easily follows from the eigenvalue problem of L on substituting Eq. 
(IA13p . Eq. (|A14p implies A\ (p,q) is made up of invariants of motion 
and T\ (p, q) is the phase of a wave in the F space that corresponds to 
the eigenvalue A. The eigenvalue problem of L is now reduced to finding 
the amplitude and the phase of waves. It is done in the main text. 
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